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COSMIC REIONIZATION AND EARLY STAR-FORMING GALAXIES: A JOINT ANALYSIS OF NEW 
CONSTRAINTS FROM PLANCK AND HUBBLE SPACE TELESCOPE 

Brant E. Robertson^ Richard S. Ellis^, Steven R. Furlanetto^ and James S. Dunlop'* 

ABSTRACT 

We discuss new constraints on the epoch of cosmic reionization and test the assumption that most 
of the ionizing photons responsible arose from high redshift star-forming galaxies. Good progress has 
been made in charting the end of reionization through spectroscopic studies of 2 ~6-8 QSOs, gamma- 
ray bursts and galaxies expected to host Lyman a emission. However, the most stringent constraints 
on its duration have come from the integrated optical depth, r, of Thomson scattering to the cosmic 
microwave background. Using the latest data on the abundance and luminosity distribution of distant 
galaxies from Hubble Space Telescope imaging, we simultaneously match the reduced value r = 0.066± 
0.012 recently reported by the Planck collaboration and the evolving neutrality of the intergalactic 
medium with a reionization history within 6 < z < 10, thereby reducing the requirement for a 
significant population of very high redshift (z ^ 10) galaxies. Our analysis strengthens the conclusion 
that star-forming galaxies dominated the reionization process and has important implications for 
upcoming 21cm experiments and searches for early galaxies with James Webb Space Telescope. 

Subject headings: galaxies: high-redshift 


1. INTRODUCTION 

Cosmic reionization represents an important era for 
assembling a coherent picture of the evolution of the 
Universe, and ambitious observational facilities are be¬ 
ing constructed to explore the most important redshift 
range 7 < z < 20. Through the Gunn-Peterson effect 
in h igh redshift QSOs and gamma ray bursts (GRBs, 
e.g., Fan et al.ll200fil Bolton et al.ll2nill : iChornock et alJ 
l2013t iMcGreer et all 1201511 and the declining visibil¬ 
ity of L yman alpha (Lyg) emission in high r e dshift 
galaxies (iStark et al.l|201 

Schenker et al.l 120121 f2(iT 


Pentericci et al. 1201 IL 12014 

Treu et al.ll201ji:lTilvi et al l 


2014f l. observations indicate that reionization ended by 
redshift z ~ 6. However, the onset and duration 
of the reionization process remain less certain. The 
most convincing constraint is provided by the integrated 
optical depth, r, of Thomson scattering to the cos¬ 
mic microwave background (GMB). The Wilkinson Mi¬ 
crowave Anisotropy Probe (WMAP) delivered a value 
T = 0.088 ± 0.014 which, in the simplest model, corre- 
spo nds to ‘instantaneou s’ reionization at Zreion — 10.5 ± 
l. KjHinshaw et ani2013ll . As a result, the WMAP result 
has been widely interpreted as i mplying that reionization 
began at z 15 or even earlier (jBromm fc Yoshid^l2011l : 
lDunlonll201,‘Hl . 

Important information on the duration of reioniza¬ 
tion can now be det ermined from the star for mation 
rate (SFR) history (jMadau fc DickinsonI 12014 here¬ 
after MD14), since early star-f orming galaxies most 
likely supply the ionizing photons (jRobertson et al.ir201fll 
1201.4) . This conclusion followed the first measures 
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of their abundance over 8 < z < 10 from Hubble 
Space Tele scope (HST) Ultra Deep Field (UDF) ob¬ 
servations (iBeckwith et afll20n6llKoekemoer et al.ll2ni,4 : 
llllingworth et al.l 1201,111 . With plausible assumptions, 
6 < z < 8 star-forming g alaxies can keep the 
Universe substantially ionized (|Robertson et al.l 120131 
iFinkelstein et al1l2Ql4) . 

However, to matc h the WMAP value of r, 
iRobertson et al.l (|201,3[ 1 also required a significant popu¬ 
lation of star-forming galaxies beyond a redshift z ~ 10. 
As a direct census of z > 10 galaxies is not currently 
possible, studies have since focused on the rate of de¬ 
cline in ab undance over 8 < z < 1 0 with mixed conclu¬ 
sions (c.f., lOesch et al.l 1201211201,4 : llshigaki et al.l [2?)T5l : 
iMcLeod et al.l 12014( 1. The requirement for a significant 
contribution of ionizing photons from z > 10 galaxies re¬ 
mains an important uncertainty whose resolution is per¬ 
ceived as a major goal for the James Webb Space Tele¬ 
scope (JWST). _ ^ _ 

The lPlanck Collaboration et al.l (j2015ll has recently re¬ 
ported a significantly lower value of the optical depth, 
r = 0.066 ± 0.012, consistent with a reduced redshift of 
instantaneous reionization, Zreion = 8.8lit'2. Here we de¬ 
termine the extent to which the Planck result reduces 
the need for significant star formation in the uncharted 
epoch at z > 10. To demonstrate this, we calculate the 
contribution of 6 < z < 10 star-forming galaxies to the 
integrated value of r, using the latest HST data. We 
then examine the residual contribution of ionizing pho¬ 
tons required from sources beyond z ~10 to match the 
new value of r from Planck, phrasing these constraints 
in terms of the likely abundance of z > 10 galaxies that 
JWST would see in a typical deep exposure. 

Th roughout we use the AB magnitude system (lOkd 
IT97I . errors represent 1 — a uncertainties, and all cos¬ 
mological calculations assume flatness and the most 
recent Planck cosmological parameters (h = 0.6774, 
Qm = 0.309, = 0.0 2230, Yp = 0.2453; 

iPlanck Collaboration et ^1201511 . 








































































2 


Robertson et al. 


2. CONTRIBUTION OF Z < 10 GALAXIES TO LATE 
REIONIZATION 

2 .1. Cosmic Star Formation History 

If Lyman continuum photons from star-forming galax¬ 
ies dominate the reionization process, an accounting of 
the evolving SFR density will provide a measure of the 
time-dependent cosmic ionization rate 

h-ion — /esc^ion/lSFRi (1) 

where /esc is the fraction of photons produced by stel¬ 
lar populations that escape to ionize the IGM, ^ion 
is the number of Lyman continuum photons per sec¬ 
ond produced per unit SFR for a typical stellar pop- 
ulation , and psfr is the cos mic SFR density. Fol¬ 
lowing iRobertson et al.l (|2013f) , we adopt a fiducial es¬ 
cape fraction of /esc = 0.2 and, motivated by the rest- 
frame U y spectral energy d istributions of z ~ 7 — 8 
galaxies (|DunloD et al.l I2013I1 . a fiducial Lyman con¬ 
tinuum photon production efficiency of logj^Q ^ion = 
53.14 [Lyc photons s~^Mq^ yr]. Somewhat larger val- 
uesof /ion may also be acceptable fe.g.. lToDDing fc Shulll 

The observed infrared and rest-frame UV luminos¬ 
ity functions (LFs) provide a means to estimate psFR- 
We use the recent compilation of IR and UV LFs pro¬ 
vided in Table 1 of MD14 and references therein to com¬ 
pute luminosity densities to a minimum luminosity 
of Lmin = O.OOIL*, where L*(z) is the characteristic 
luminosity of each relevant LF parameterization (e.g., 
Schechter or broken power law modelsj^- We supple¬ 
ment the MD14 compilation by including psFR values 
computed from the LF determinations at z 8 by 
iSchenker et al.l (|2013ll . at z ^ 7 — 8 bv iMcLure et ^ 
(I2013|l^ and estimates at z ^ 10 by lOesch et al.l (I20l3 l 
and iBouwens et all (|2014ll . We inclu de new HST Fron- 
tier Fields LF cons Uaints at z ^ 7 bv lAtek et al.l (120141 1 
and at z ^ 9 by iMcLeod et a,l.l (1201411 , incorporatin g 
cosmic variance estimates from IRobertson et alJ (|2014ll . 
We also updat e d the MD14 estimates derived from the 
IBouwens et aLI (1201211 LFs at z ^ 3 — 8 with newer mea- 
surements bv bouwens et al.l (|2014ll . All data were con¬ 
verted to the adopted Planck cosmology. 

We adopted the conversion psFR = kpl supplied 
by MD14 for IR and UV luminosity densities, i.e. 
Kir = 1.73 X 10“^° Mq yr“^ Lq^ and kuv = 2.5 x 
10“ Mq yr Lq^ respectively, as well as their redshift- 
dependent dust corrections and a Salpeter initial mass 
function. Uncertainties on psFR are computed using 
faint-end slope uncertainties where available, and other¬ 
wise we increased the uncertainties reported by MD14 
by the ratio of the luminosity densities integrated to 
L = 0.03L* and L = O.OOIL*. The data points in Fig¬ 
ure [T] show the updated SFR densities and uncertainties 
determined from the IR (dark red) and UV (blue) LFs, 
each extrapolated to Lmin = O.OOIL*. 

Since we are interested in the reionization history both 
up to and beyond the limit of the current observational 
data, we adopt the four-parameter fitting function from 

® We ado pt this limit sinc e it co rresponds to Mmax ~ —13 at 
2 ~ 7, which lRobertson et all II2013I 1 found was required to reionize 
the Universe by 2 ~ 6. It corresponds to Mmax = M* + 7.5. 



Fig. 1.— Star formation rate density psFB, wi t h red shift. Shown 
are the SFR densities from IMadau fc DickinsonI Il2014l ~) determined 
from infrared (dark red points) and ultraviolet (blue points) lumi¬ 
nosity densities, updated for recent results and extrapolated to a 
minimum luminosity Lmin = O.OOIL*. A parameterized model for 
the evolving SFR density (Equation [2]l is fit to the data under the 
constraint that the Thomson optical depth t to electron scatter¬ 
ing measured by Planck is reproduced. The maximum likelihood 
model (white line) and 68% credibility interval on pspR (red re¬ 
gion) are shown. A consistent SFR density history is found even 
if the Planck t constraint is ignored (dotted black line). These 
inferences can be compared with a model forced to reproduce the 
previous WMAP t (orange region), which requires a much larger 
PSFR 3-t redshifts 2 > 5. 

MD14 to model psfr(z). 


Psfr(2:) 


_ (1 + 

^^1 + [(1 -b ^)/Cp]‘^P 


( 2 ) 


and perform a maximum likelihood (ML) determination 
of the par ameter values us ing Bayesian methods (i.e.. 
Multinest: iFeroz et al.l [200^ assuming Gaussian errors. 
If we fit to the data and uncertainties reported by MDI4, 
we recover similar ML values for the parameters of Equa¬ 
tion [21 The range of credible SFR histories can then be 
computed from the marginalized likelihood of psFR by 
integrating over the full model parameter likelihoods. 


2.2. Thomson Optical Depth 

If photons from star forming galaxies drive the reion¬ 
ization process, measures of the Thomson optical depth 
inferred from the CMB place additional constraints on 
PSFR- The Thomson optical depth is given by 

t{z) = c{nn)(TT [ ( 3 ) 

where c is the speed of light. The comoving hydrogen 
density (nn) = Xp^ltPc involves the hydrogen mass frac¬ 
tion Ap, the baryon density 17;,, and the critical density 
Pc- The Thomson scattering cross section is ut- The 
number of fre e electrons per hydrogen nucleus i s calcu¬ 
lated following iKuhlen fc Faucher-Gigubi^ (|2012D assum¬ 
ing doubly ionized helium at z < 4. 
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Fig. 2.— Thomson optical depth to electron scattering r, in¬ 
tegrated over redshift. Shown is the Planck constraint r = 
0.066 ib 0.012 (gray area), along with the marginalized 68% credi¬ 
bility interval (red region) computed from the SFR histories psFR 
shown in Figure [T] T he corresponding inferences of r{z) from 
IRobertson et afl 11^131) (dark blue region), a model forced to re¬ 
produce the 9-year WMAP r constraints (orange region), and a 
model with psfr truncated at 2 > 8 (light blue region) following 
lOesch et al.l l|2014l ^ are shown for comparison. 

The IGM ionized fraction is computed by 

evolving the differential equation 


Qhii 


hion Chii 

(uh) tree 


(4) 


where the IGM recombination time 


tree = [CH„aB(T)(l + Vp/dVp) (uh) (1 + (5) 

is calculated by evaluating the case B recombination co¬ 
efficient ae at an IGM temperature T = 20, OOOK and 
a clumping fract ion Chu = 3 fe.g.. iPawlik et all 120091 : 
iShull et al.ll2012|l . We incorporate the Planck Thomson 
optical depth constraints (r = 0.066 ±0.012, treated as 
a Gaussian) by computing the reionization history for 
every value of the psFR model parameters, evaluating 
Equation O and then calculating the likelihood of the 
model parameters given the SFR history data and the 
marginalized Thomson optical depth. 

Figure [T] shows the ML and 68% credibility interval 
(red region) on /0sfr(-2) given the psFR constraints and 
the newly-reported Planck Thomson optical depth. We 
find the parameters of Equation [5] to be Op = 0.01376 ± 
0.001 Mq yr Mpc"^ bp = 3.26 ± 0.21, Cp = 2.59 ± 0.14, 
and dp = 5.68 ± 0.19. Without the Thomson optical 
depth constraint, the values change by less than 1%. 
These inferences can be compared with a SFR history 
(Figure [1] orange region) forced to match the previous 
WMAP measurement (r = 0.088 ± 0.014) by upweight¬ 
ing the contribution of the derived r value relative to the 
PSFR data. The model’s ML parameters (op = 0.01306, 
bp = 3.66, Cp = 2.28, and dp = 5.29) lie well outside 
the range of models that reproduce jointly psfr(-z) and 


the Planck r. Fitting to only data at z > 3 or only in¬ 
dependent data points at z > 6 changes our credibility 
intervals by ~ 25%. 

We can now address the important question of the 
redshift-dependent contribution of galaxies to the Planck 
r = 0.066 ± 0.012 in Figure [21 The red region shows 
a history which is consistent with the SFR densities 
shown in Figure |T] given our simple assumptions for the 
escape fraction /esc, early stellar populations, and the 
dumpiness of the IGM. Importantly, the reduction in r 
by Planck (compared to WMAP) largely eliminates the 
tension between pc;fb(z:) and r that was discu ssed by 
many authors, including IRobertson et al.l (1201311 . That 
a SFR history consistent with the psfr(-z) data easily 
reproduces the Planck r strengthens the conclusions of 
IRobertson et al.l (|2013fl that the bulk of the ionizing pho¬ 
tons emerged from galaxies. Figure [2] shows that the ob¬ 
served galaxy population at z < 10 can easily reach the 
68% credibility intervals of r with plausible assumptions 
about /esc and Amin- As a consequence, the reduced r 
eliminates the need for very high-redshift (z 10) star 
formation (see section 3 below). We note the dust cor¬ 
rection used in computing psFR at z ^ 6 permits an 
equivalently lower /esc without significant change in the 
derived r. We note that to reach r > 0.08 given the 
Psfr(- 2) constraints requires /esc ^ 0.3 or Chu ^ 1- 

Figure [2] also shows r(z) computed with the 9—year 
WMAP T marginalized likelihood as a constraint on the 
high- redshift SFR density (blue region; IRobertson et al.l 
1201311 . which favored a relatively low r ~ 0.07. If, in¬ 
stead, the SFR density rapidly declines as psfr oc (1 ± 
z)~^°' ^ beyond z ~ 8 as suggested by, e.g., lOesch et al.l 
(1201411 . the Planck T is not reached (light blue region). 
Lastly, if we force the model to reproduce the best-fit 
WMAP T (orange region), the increased ionization at 
high redshifts requires a dramatic increase in the z > 7.5 
SFR (see Figure |T]) and poses difficulties in matching 
other data on the IGM ionization state, as we discuss 
next. 


2.3. Ionization History 

Similarly, we can update our understanding of the 
evolving ionization fraction (3Hn('2) computed during 
the integration of Equation |4l Valuable observational 
progress in this area made in recent years exploits 
the fraction of st ar forming galaxie s showing Lyman- 
a emission (e.g., iStark et al.( 1201(111 n ow extended to 
z ^ 7 — 8 from iTreu et ^ j20l5l . iPentericc et al.l 
(|2014ll and lSchenker et ahT i 201411 . the Lyman-a damping 
wing absorption const raints from GRB host galaxies by 
iChornock et al.l (|2013[ 1. and the number of dark pixels 
in Lyman-g forest ob servations of background quasars 
(iMcGreer et al.l [20151 1. While most of these results re¬ 
quire mo del-dependent inferences to relate observables 
to Qhii, they collectively give strong support for reion¬ 
ization ending rapidly near z ~6. 

Figure |3| shows these constraints, along with the in¬ 
ferred 68% credibility interval (red region; ML model 
shown in white) on the marginalized distribution of the 
neutral fraction 1 — Qhh from the SFR histories shown 
in Figure |T] and the Planck constraints on r. Although 
our model did not use these observations to constrain the 
computed reionization history, we nonetheless find good 
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ML + 68% Credibility Interval 
Robertson et al. 2013 
Forced Match to WMAP r 
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Fig. 3.— Measures of the neutrality 1 — Qhh inter- 

galactic medium as a function of r edshift. Shown are th e ob¬ 
servational constraints compiled by IRobertson et all 1)20131 ') , up¬ 
dated to include recent IGM neutrality est imates from the ob¬ 
served fraction of Lym an-a emitting galaxies JSchenker et al.l2014l : 
IPenteric cie^aj|20^), constrai nts from the Lyman-a of GRB host 
g:alaxies~ }^hornock et ^120131) , and inferences from d ark pixels in 
Lyman-a forest measurements JMcGreer et al.ll2Q15ll . The evolv¬ 
ing IGM neutral fraction computed by the model is also shown (red 
region is the 68% credibility interval, white line is the ML model). 
While these data are not used to constrain the models, they are 
nonetheless remarkably consistent. The bottom panel shows the 
IGM neutral fraction near the end of the reionization epoch, where 
the presented model fails to capture the complexity of the reion¬ 
ization process. For ref erence we also show the corresponding in¬ 
ferences calculated from [Robertson et al.l II2013I1 (blue region) and 
a model forced to reproduce the WMAP r (orange region). 

agreement0. 

Figure a ls o sh ows the earlier model of 
IRobertson et al.l (|2013ll (blue region) which com¬ 
pletes reionization at slightly lower redshift and displays 
a more prolonged ionization history. This model was in 
some tension with the WMAP r (Figure [2]). If we force 
the model to reproduce the WMAP t (orange region), 
reionization ends by z ^ 7.5, which is quite inconsistent 
with several observations that indicate neutral gas 
within IGM over the range 6 < z < 8 (Figure |31). 

3. CONSTRAINTS ON THE CONTRIBUTION OF Z > 10 
GALAXIES TO EARLY REIONIZATION 

By using the parameterized model of MD14 to fit the 
cosmic SFR histories, and applying a simple analytical 
model of the reionization process, we have demonstrated 
that SFR histories consistent with the observed psfr{z) 
integrated to = O.OOIL* reproduce the observed 

Planck r while simultaneously matching measures of the 
IGM neutral fraction at redshifts 6 < z < 8. As Fig¬ 
ure [T] makes apparent, the parameterized model extends 

® The model does not fare well in comparison to Lyman-a forest 
measurements when Qhu ~ 1 because of our simplified treatment 
of th e ionization process (see the discussion in IRobertson et aEl 
120131 


Fig. 4. — Correspondence between the Thomson optical depth, 
the equivalent instantaneous reionization redshift ^reioni ^Tid the 
average SFR density psFR redshift z: > 10. Shown are samples 
(points) from the likelihood function of the pSFR model parameters 
resulting in the 68% credibility interval on r from Figure [J] color 
coded by the value of ^reion- The samples follow a tight, nearly 
linear correlation (dashed line), demonstrating that in this model r 
is a proxy for the high-redshift psFR- We also indicate the number 
of z: > 10 galaxies with ruAB < 29.5 per arcmin”^ (right axis), 
assuming the LF shape does not evolve above z > 10. 

the inferred SFR history to z > 10, beyond the reach 
of current observations. Correspondingly, these galaxies 
supply a non-zero rate of ionizing photons that enable the 
Thomson optical depth to slowly increase beyond z ^ 10 
(see Figure [2|). We can therefore ask whether a connec¬ 
tion exists between psfr(,z > 10 ) and the observed value 
of T under the assumption that star forming galaxies con¬ 
trol the reionization process. 

Figure 2] shows samples from the likelihood function of 
our model parameters given the psfr)^:) and r empirical 
constraints that indicate the mean SFR density (psfr) 
(averaged over 10 < z < 15) as a function of the total 
Thomson optical depth r. The properties (psfr) and r 
are tightly related, such that the linear fit 

(psfr) ~ 0.344(r - 0.06) -b 0.00625 [Mq yr"^ Mpc“^] 

( 6 ) 

provides a good description of their connection (dashed 
line). For reference, the likelihood samples shown in Fig- 
urelHindicate the corresponding redshift of instantaneous 
reionization Zreion via a color coding. 

Given that the SFR density is supplied by galaxies 
that are luminous in their rest-frame UV, we can also 
connect the observed r to the abundance of star form¬ 
ing galaxies at z > 10. This quantity holds great in¬ 
terest for future studies with James Webb Space Tele¬ 
scope, as the potential discovery and verification of dis¬ 
tant galaxies beyond z > 10 has provided a prime moti¬ 
vation for the observatory. The 5-cr sensitivity of JWST 
at 2 ^m in a t = 10^ s exposure is toab ~ 29.50 

See http://www.stsci.edu/jwst/instruments/nircam/sensitivity/table 
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At z ~ 10, this sensitivity corresponds to a UV ab¬ 
solute magnitude of Muv ~ —18. Extrapolating the 
SFR density to 2 > 10 and using the shape of the LF 
at z > 9, we estimate that N ~ 0.5 arcmin”^ galax¬ 
ies at z > 10 will be present at apparent magnitudes of 
TOAB < 29.5 at A = 2 /rm. Deep observations with JWST 
over ~ 10 arcmin ^ may therefore find 5 candidates at 
z > 10 (see also iBehroozi fc 81151201511 . Returning to 
Figure [H we can see the impact of the reduced value of 
r by comparing the Planck and WMAP curves beyond 
z ~ 10. 

4. DISCUSSION 

The lower value of the optical depth t of Thomson 
scattering reported by the Planck consortium (2015) 
strengthens the likelihood that early star-forming galax¬ 
ies dominated the reionization process, as our model can 
simultaneously match the observed SFR history (Figure 
[1]) over 6 < z < 10, the integrated value of r (Figure [2]), 
and recent constraints on the IGM neutral fraction over 
z~6-8 (Figure®. 

A state-of-th e- art reionization analysis by 
iChoudhurv et al.l (1201411 used the distribution of 
Lyman a equivalent widths, the IGM photoionization 
rate, and the mean free path of ionizing photons, to also 
conclude that reionization likely comp leted at z ^ 6, 
with a corresponding r ss 0.07(see also iRobertson et'all 
[Mil . With Planck now favoring r ss 0.066 and 
informed by a full accounting of available constraints on 
the SFR history, we have reached similar conclusions 
using different empirical inputs. 

Our modeling makes some simplifying assumptions, 
adopting a constant escape fraction /esc = 0.2, IGM 
clumping factor C Ri 3, and Lyman continuum pro- 
duction efhciencY for early stellar populations. In 
IRobertson et al.l (|20I3ll we examined these assumptions 
carefully and tested more complex models, e.g. with 
evolving escape fraction require d to match the IGM pho - 
toionization rates at z < 6 (e.g.. iBecker fc Boltonll20lll . 
These assumptions influence the computation of r and 
Qhii but do not affect the inferred SFR history in Figure 
dl Our conclusion that z < 10 galaxies can account for 
the Planck r relies on extrapolating LFs below observed 
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limits and a higher escape fraction than at lower redshift. 
If galaxies are less efficient ionizers, mor e z > 10 star for¬ 
matio n would be permitted. However, IRobertson et all 
(l20I3fl already demonstrated such an ionizing efficiency 
is required to maintain a highly ionized IGM at z ~ 7 
(Figure®. 

The “excess” value of r above that provided by galax¬ 
ies at z < 10 measures psFR at z > 10. Equation [6] and 
the Planck I-tr upper limit on r provide an upper limit 
of Psfr(2 > 10) < 0.013 Mq yr“^ Mpc“^. This provides 
the first empirical limit on models that increase the ion¬ 
izing efficiency during this epoch e.g. with massive Pop- 
ulation III stars and sta r formation in mini-halos (see 
iLoeb fc Furlanett^ 120131 for an overview of such mod¬ 
els). Our results suggest such models cannot dramati¬ 
cally change the star formation efficiency at the earliest 
times. 

Reionization proceeds relatively quickly as the ionized 
fraction evolves from Qhh = 0-2 to Qhh = 0-9 in only 
400 Myr of cosmic history over 6 ;< z < 9. This du¬ 
ration is consistent with recent up per limits on the k i- 
netic Sunyaev-Zel’dovich effect (e.g.. lGeorge et al1l2015[l . 
Our results offer extra hope for efforts to make red- 
shifted 21-cm measurements of neutral hydrogen in the 
IGM, as the ex perimental foregrounds ar e weakest at low 
redshifts (e.g.. Ivan Haarlem et ahi l2013t iBowman et al.l 
120131 : iPober et al.ll20lJl . Such experiments are essential 
for testing key assumptions in our analysis (like /esc and 
Amin) by observing the reionization process directly. The 
apparent lateness of reionization suggests that next gen¬ 
eration experiments, which hope to reach z 20, can 
probe even earlier phases of galaxy formation. 
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